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■ The simple algorithm for the simulation and visualization of non relativistic 

quantum dynamics is proposed that is based on a collective behavior of classical 
Q_i! particles. Any quantum particle is represented as the swarm of its classical sam- 

ples which interact by simple rules including emission and absorption of samples of 
tied photons. The quantum dynamics results from the collective behavior of such 
^ ' a swarm where the eigenstates are treated as the equilibrium states relatively to 

emission-absorption of photons. The entanglement is treated as a correlation be- 
tween samples of the different swarms that is stored in the space-time part of the 
^ , model inaccessible for a user. The amplitude is always grained. The Coulomb field 

^ ' between quantum particles is simulated, analogously to free flow of quantum pack- 

age, by the point wise interaction between its samples and scalar photon samples 
which propagate by diffusion. It gives square root speedup in comparison to each- 
with-each method. This method obviously includes decoherence and admits the 
natural generalization on the QED of many particles with the linear computational 
cost. 
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1 Introduction 



Many particle problems represent the natural area for the applications of quantum theory, 
and just here it meets the unexpected and principal difficulty At first glance this difficulty 
seems not fatal because it is connected with the complexity of the description of many 
body problems on classical computers, e.g., with the applications of quantum physics. But 
this viewpoint is wrong. In the more detailed consideration the more and more depths of 
this difficulty become apparent. We can now assert with certainty that this difficulty is so 
serious that it factually establishes the limitation of the applicability of quantum theory 
itself, so that the adequate description of many body systems requires the substantial 
revision of its basement: namely, of its mathematical apparatus. The difficulty is that 
the dimensionality of the space of states of many body systems grows as the exponential 
when the number of its particles increases. The natural solution would be the creation of 
a quantum computer for the simulation of many body quantum systems. 1 . The idea of 
the creation of quantum computer obviously belongs to the great ideas of physics. The 
valuable part of physicists (including the author of this paper) is somehow involved in 
this work. The quantum computer is the first principle prediction of quantum theory 
that follows immediately from the mathematical formalism of tensor products of Hilbert 
spaces. Hence, the creation of such a device must be treated not as the technological 
task but as the decisive experiment on the checking of quantum theory in the many 
particle case. It worth noting that quantum theory has been never verified in many 
particle case; all the results experimentally confirmed at the moment are obtained by 
the different reduction to the one particle case 2 The factual implementation of quantum 
theory in the many body area, e.g., the practical, computational consideration of the wave 
functions of the form ip(ri,r 2 , . . . ,r n ) for large n cannot be realized neither analytically 
nor numerically due to the phenomenon of quantum speedup of classical computations 
that was discovered in ninetieths. The sense of this phenomenon is as follows. Accordingly 
to Hilbert formalism, the states of system with n particles belong to the tensor product 
Tii (^) 7~^2 ® • • • ® 7~Ln of spaces Tij for one particle states (for simplicity we assume that 
all particles are distinguishable; otherwise we should consider Fock space of occupation 
numbers with bosons or fermions symmetry instead of Hilbert space that complicates the 
representations without changing the essence). And conversely, any normalized vector in 

1 The idea of quantum computer was in the air and has been expressed almost simultaneously by few 
people: Benioff, Manin, later - Deutsch, but Feynman was the first who formulated it clearly as applied 
to the simulation of quantum many body problems, he also reformulated quantum theory on the intuitive 
language of paths (that corresponds for our aim better than the traditional Shroedinger form). Feynman, 
of course, did not state the question about the limitation of applicability of Hilbert spaces apparatus for 
that it was not necessary experimental data at that time. 

2 The single exclusion is the attempts of the direct computation of the energy of bipartite systems via 
wave function ip{r\,r2), for example, for helium atom (see below). The attempts to solve many body 
problems analytically always ended by ambiguous results: for example, the systematic consideration of 
evolution of two photons and one electron in QED leads to the divergence of rows for the amplitude. 
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this space represents some physically realizable state of n particle quantum system. The 
dimensionality of this space equals the product of dimensionalities of all Hj, e.g., it grows 
rapidly with the growth of n. It means that there exist such states that amplitudes of 
some basic states in it are very small. It is impossible to check the presence of such small 
amplitudes by the straightforward measurements because the corresponding probabilities 
are negligible. 

But there exist the special forms of many body system evolutions when these small 
amplitudes constructively add in the huge quantities and give the observable value in 
so short time that no classical method of computations (in particular, no analytical cal- 
culations) can keep pace with this process. There are the quantum computations. It 
is important that there exists the exact criterion separating a quantum computer from 
classical ones: that is the speed of solution of some computational tasks (for example, the 
search problem). We thus have the project of decisive experiment checking the truth of 
quantum theory in the many body area. In the present time all the attempts to build at 
least a few qubit reliable quantum computers (let along a scalable one) failed. Of course, 
one cannot yet conclude that the quantum theory cannot be applied to the many particle 
tasks. But it makes the building of alternative formalism for many body problems actual. 
Such formalism must embrace all known quantum physics including entangled states but a 
scalable quantum computer. The single possible mathematical formalism for such extrap- 
olation of quantum, theory to the many particle case is the classical algorithm theory. The 
simplest idea is to introduce the grain or the quantum of amplitude (the minimal nonzero 
module of amplitude) for quantum states and to consider the corresponding reduction of 
quantum theory, which was proposed in the work |3] 3 . It turns that this simple trick gives 
the unified description of quantum evolutions without division to the unitary evolutions 
and the measurements (including decoherence) . The Born rule for quantum probability 
follows from the hypothesis of amplitude quantum immediately. But this method is not 
yet ready for practical implementations. It is based on the notion of many-body wave 
function and for the computer realization it requires the immediate using of the amplitude 
quantum which real value can turn too small for this. It causes over expenditure of the 
computational resources and thus allows solving mainly one particle problems. 

In this work we consider the version of algorithmic approach (see. |3]), based on the 
collective behavior of classical many body systems. This version can be named swarm 
method, because it uses the representation of a quantum particle as a swarm of its samples, 
each of which is a classical particle and can interact with others by simple rules. Any 
set of identical quantum particles is represented by the corresponding joint swarm. This 
method admits the easy passage from the classical description of particle to the quantum 
and vice versa, it is scalable and convenient for the creating of computer programs. In 
this version only EPR type of entanglement can be adequately represented that is the 

3 It makes our model dependent on the choice of basis. This cost seems not too large because algebraic 
universality is a feature of Hilbert spaces formalism and it does not follow straightforwardly from the 
experiments. 
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essential reduction of Hilbert quantum formalism. In the framework of swarm method the 
fundamental interaction between charged particles and photons can be interpreted, that 
means that this method can embrace QED, and probably some other types of interactions. 

2 General features of algorithmic approach: user and 
administrative parts of model 

We make some overview of algorithmic approach. The main tool here is classical algo- 
rithms, not the analytical technique. The final result is represented by the algorithm 
simulating the dynamics of system at hand, not by the value of some parameters. From 
the theory of algorithms it is known that the work of algorithm cannot be predicted be- 
forehand. The only way to learn how it works is to launch it and wait when it returns 
the result. The single universal form of the evaluation of results is thus the observation 
of video film which is the visualization of the work of the simulating algorithm. The 
usual physical methods like finding Hamiltonians, eigenvalues and eigenvectors etc., will 
obtain the status of auxiliary instruments, checkpoints or heuristics in building of the 
final dynamical picture of evolution. This results in the division of the model to the user 
and the administrative parts. The first one is accessible for user (it inserts the initial data 
and shows the video film). The second part is inaccessible for a user; it prepares the film 
through the simulating algorithm. 

The physical sense can be acquired for objects containing in the user part only, which 
can be observe immediately. As for the administrative part, its details have no physical 
sense. The administrative part is needed for the film preparation only and it cannot be 
associated with some object in the real world. In particular, the time spent to the film 
preparation is not the physical time 4 . The necessity of the administrative part follows 
from the existence of quantum non locality, e.g., entanglement of spatially separated 
particles. The simulation of such systems cannot be fulfilled in the real time mode, the 
film must be prepared beforehand to be able to work at the both parts of a spatially 
separated system in an entangled state. Here it would be naturally to provide by the 
free will only users, but not the administrative part of the model 5 . For example, if 
we need a source of random numbers we can include it to the administrative part as 
some program generating pseudo random numbers. The seeming absolute character of 
quantum randomness will thus result from that we do not know the parameters of this 
generator which lies in the administrative part. We also should not be excited by the 
possibility to learn these parameters by some non authorized way and then use it for 
our technologies. The right simulating algorithm produces the film which is independent 
from the individual parameters of randomness generator. The radical viewpoint is that 
the randomness generators are needed for the debugging of the true model embracing 

4 Some authors use the term hidden time of quantum theory, see, for example, 
5 I do not intend to go into the details of the philosophic character. 
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the huge number of particles, for example, in biology. We will use the generators of 
randomness because our immediate aim is the physics of many-body systems where the 
randomness is the fundamental element of evolution. 

The bipartite structure of the model would make no sense without the imposing of 
strict limitations to the abilities of the administrative part. This limitation is only one our 
model must represent the classical algorithm realizing effective (polynomial complexity) 
algorithms. If we consider the whole Nature as such a model, we would treat ourselves 
as its users and the administrative part as the internal structure of the world which im- 
mediate access is forbidden for us, for example, by some password. Despite of the closed 
character we can make one conclusion about its abilities: it cannot realize fast quantum 
algorithms. This feature of the model allows the contra version and it can be experi- 
mentally checked. Such experiment is the building of a scalable (10000 or more qubits) 
quantum computer. The creation of such device would mean the crash of algorithmic 
approach. But no failure in quantum computer building can serve as the proof of cor- 
rectness of algorithmic approach. The correctness of complex films can be the unique 
criterion here. 

The main aim in the algorithmic approach is the creation of films reflecting the real 
evolutions of many-body quantum systems. Just this aim we keep in mind describing our 
model. As for the correspondence with the standard quantum physics, such a correspon- 
dence cannot be complete, because we at first exclude the description of n particle system 
via wave function of the form \E , (r 1 ,r 2 , . . . ,r n ). This description is possible in our model 
only for the case of the entangled states of EPR type, for example, one photon and one 
particle, two entangled photons, or many particles which behave as the single particle, 
etc. In the reality this assumption is not very restricting; we can find substantially exact 
corresponding with QED that is the subject of future investigations. We here take up 
only non relativistic approximation. Nevertheless, it will be clear how to find the rules 
for photon emission- absorption which agree with QED. 

3 Swarm model of quantum dynamics 

Swarm model (or model of collective behavior) consists of that each quantum particle is 
represented as the set (swarm) of classical particles of the form 

s = {s 1 ,s 2 ,...,s k }, (1) 

where each particle Sj is called a sample of the considered quantum particle and has the 
definite space-time coordinates t(sj),x(sj) and some parameters needed for the simulation 
of dynamics. For example, we can introduce for each sample the pointers to the other 
samples of this quantum particle or the other quantum particles. The pointers associated 
with Sj point to the samples which influence is most essential for Sj. In the swarm 
representation of charged particles space-time coordinates are separated because these 
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particles are classical and their impulses are always connected with the coordinates as 
P = m l±[- ^ ne c l° se coordinates in QED is the good criterion of the essentiality of 
interaction, and consequently we can manage without the pointers at all. The simple 
swarm representation of n charged particles has thus the form of list s — (s 1 , s 2 , . . . , s n ), 
which components s k are swarms representing the separate particles. This simple swarm 
representation is valid for the description of the state of this system in any time instant. 
In the language of standard Hilbert formalism it means that this state will always has the 
form (^) v I / 2 <S) • • • ® ^n, e.g., will be not entangled. This is the fundamental difference 
between swarm method and Hilbert formalism, because the most essential property of 
the last is the existence of entangled states. It worth noting that we could simulate the 
standard formalism via pointers from the samples of charged particles to the other samples 
of this sort. 

If we introduce the operator of free flight for the swarm representation and some simple 
rules like creation-annihilation of samples, the simple swarm representation is valid for 
finding of the ground states of a system with n electrons (the anti symmetry of wave 
function can be reflected in the framework of this representation). But such operator 
cannot simulate the dynamics of quantum systems even in case of one particle in potential 
(excluding the case of free particle flow which is simulated by the diffusion of its samples). 

To simulate the unitary evolution induced by Shroedinger equation we need to intro- 
duce the samples of the different sort. These samples must 1) move very fast, 2) they must 
be emitted by the samples of charged particles, and they must transform to the samples 
of charged particles, and 3) they must store the memory about what sample of charged 
particle they are emitted from; e.g., there must exist the pointers connecting them with 
samples of charged particles. In view of these properties we call these new samples the 
samples of tied photons. We can introduce vector photons which direction of propagation 
is orthogonal to the vector E of electric field, and the scalar photons moving along E. 
At last we can define the conditions of transformation of samples of tied photons to the 
samples of free photons that is needed for representation of QED. 

The pointer of sample a of particle or photon sample 7 is the natural number p(a) 
or ^(7), denoting the number of the other sample - of photon or particle correspondingly. 
The pointers must satisfy the reciprocal condition: pi^k) = m if f p(a m ) = k. 

The swarm representation of n particle system in electromagnetic field has the form 
of list of vectors 

Spar = {s\s 2 ,...,s n }, 7 = {71, 72, • • • , 7z}, <7 = {#1, 02, • • • , 0d}, (2) 

where s par is the vector consisting of samples of charged particles with pointers particle- 
photon to the samples of photons which form the vectors 7, 9 of the tied and free photon 
samples. The evolution of this representation (J2J) in time is determined by the sequential 
operators of the following types: creation of samples of tied or free photons by the samples 
of charged particles (photon emission), creation or annihilation of the samples of charged 
particles by the photon samples (absorption) - here the sample of tied photon disappears 
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and the sample of free photon remains, and also creation and annihilation of the samples 
of charged particles and photons. The interaction of QED can be simulated by this type 
of interactions. 

From the viewpoint of Hilbert formalism a state of system of n quantum particles we 
deal with in swarm method will always have the form 

|tf) = |*l>(g)l*2>(g)...|tfn>, (3) 

where G TLparticies *S)7~Lj, photons is not entangled state of j-th particle and photons of 
the form 

^ j ^ j, par (^^) ^3 phot (4) 

phot is the state of some subset Phj of all photon samples so that the sets Phj, j = 
1,2, ... ,s form the division of the set Ph of all photon samples. No entanglement thus 
arises in the passage to the standard language of Hilbert spaces. 

Here the operator of time evolution for such system has the following form: 

^* a ,*3,...,*« ^ ; (g) . . . (g) I/* 1 '**--*"- 1 , (5) 

where j+i>->*« j s ^he superposition of a) unitary operator on j-th particle 

which depends on the states of all other particles and b) creation and annihilation of 
photon samples leading to the change of the sets Phj. We thus have only not entangled 
states of particles and photons but an evolution of such states cannot be represented as 
sequence of one particle operators. The evolution induced by the operator (jSJ), cannot 
be expressed in form of quantum gate array because the application of any entangling 
operator gives the entangled state. The method of collective behavior for many particles 
can be thus expressed in terms of tensor products, e.g., here we have the radical difference 
from the many body Hilbert formalism. Operators of the form (jSj) arise in mean field 
approximation that is applied in the Hartree-Fock method. This method and the diffusion 
Monte Carlo method give good approximations of the experimental values of energy for 
many electron systems (atoms and molecules), with the accuracy within few percents. 
The most efficient methods of many body computations thus can be well approximated 
by thhe swarm method. 

The entangled states of the form: 

j 

have satisfactory swarm representation, where for all k = l,2,...,n . . . is 

orthonormal system of functions in the space H,k of states of particle or photon k. Such 
states generalize the known Schmidt decomposition for bipartite states but only in two 
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particle case (n = 2) any state can be reduced to them. Swarm representation of the state 
© has the form (J2J) where s k is the swarm representing the state 

5>i)i*j>, ( ? ) 

3 

where pointers are arranged so that they connect samples corresponding to the basic one 
particle functions \*f>j) with the same lower index j. 

Swarm approach thus principally reduces the possible types of entanglement which 
exist in tensor products of Hilbert spaces. 

4 Swarm representation of many-body wave func- 
tions as alternative to Hilbert formalism 

The principal difference between the method of collective behavior and Hilbert formalism 
is that in the first method the four dimension space-time is the single space of states and 
its dimensionality does not grow with the number of particles. It results from that the set 
of quantum particles is represented by the set of the corresponding swarms which exist 
in the same space-time. The versions of swarm approach are the known diffusion Monte 
Carlo method and its modifications like Semenihin method for finding of excited states 

(El)- 

Our aim here is the representation of the quantum dynamics of a system consisting of 
many charged particles interacting through photons accordingly to QED. 

4.1 Representation of one particle wave function through col- 
lective behavior with tied photons 

We give more universal representation of quantum evolution which can be then generalized 
to the case of many particles. Shroedinger equation is twice more complex than the 
diffusion transform because the wave function *&(x, t) - is complex- valued. It carries 
twice more information about the state of swarm than the density p from the diffusion 
equation, because besides coordinates of samples we must know their speeds. We are 
going to describe one particle quantum dynamics in terms of swarm method. The simple 
considerations show that it is impossible if we are limited by only samples of one particle 
even if we will divide them to the different types; arguments for that see below. Hence, 
we must introduce the samples of some other auxiliary particle (tied photon), different 
from the initial. The role of this auxiliary particle is to ensure the quantum dynamics 
of the initial one. Here the stationary functions \<f>j) of Hamiltonian correspond to the 
equilibrium states sf at of the swarm s, that means the steadiness of necessary storage of 
the tied photon samples, whereas in the passage from one stationary state to another we 
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need to change this storage by increasing or decreasing it. It makes reasonable to call this 
new particle tied photon, despite of that such characteristics of photons as polarization 
will not be yet included to the description of simple quantum dynamics. In this section 
the association of new particle with electrodynamics photon will have no consequences 
because we will consider not a real particle but its samples only. 
Shroedinger equation for the wave function *&(x, t) has the form 

ii = -A* + V*. (8) 
We represent the wave function in the form 

V(x,t) = %? r (x,t) + i¥(x,t), 

where ty r (x, t), ^ l (x,t) are its real and imaginary parts. The equation (jHJ) can be then 
written as the system of two equations 

& = -A¥ + V¥, , , 



We establish the necessary connection between Shroedinger and diffusion equations. 
The diffusion equation has the form 



p\ii = div(p grad u) — qu + F, (10) 

where u(x, t) is the density of particles (samples), p%,p, q, F parameters depending on x, t 
and designing the density of environment, diffusion coefficient, coefficient of absorption 
and intensity of source of samples correspondingly. The positive coefficient of absorption 
means that the samples are annihilated in this point with the intensity \q\. We assume 
that F — 0, and p\,p has unit value, e.g., we will give the reasonable sense only to 
the coefficient of absorption q, which will be proportional to the potential energy in 
Shroedinger equation. The equation (fTUI) follows from Nernst law (see. (j7]) for the 
finding of the string of samples through the element of surface dS: 

fin 

dQ = -p—dS. (11) 
on 

To reduce the equation (jHJ) to some version of the diffusion equation, we consider the 
swarm of samples of our quantum particle. We divide these samples to two types: real (r) 
and imaginary (i), and in each of these type in turn divide to two subtypes: positive (+) 
and negative (-). We thus obtain the division of all samples to four types which members 
will be designated as follows: cnp r , o;p\ ctj' r , ctj'\ where j denotes the number of a 
sample. For the description of stationary states only one type of samples would suffice, 
because stationary states are determined by the density. For the description of dynamics 
we really need two types of samples: real and imaginary, but it is more convenient to 
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have four types. We divide the configuration space to the small cubes and let D(x,t) 
denote the cube containing the point (x,t). The total number of samples of the swarm 
s of a same type containing in some cube D(x,t) will be denoted as s cr,v (x,t), where 
a G {+, — }, r\ G {r, i}. We assume that the speeds of samples are distributed uniformly 
independently of the type. We intend to represent the wave function evolution as the 
chain of sequential diffusion we must get rid of signs in equations, and for that we will use 
the defined types of samples where the swarm approximation of wave function is always 
found by the formula 



*{x,t) s 



(12) 



This equation does not determined uniquely the division to the positive and negative 
parts, but only within the addition of a constant to these both parts. We then have 
within the normalization the following approximate equations: 



(x,t) 
{x,t) 



X, t) 



s + ' l (x, t) — S~ 



r (x,t), 

Kx,t), 



(13) 



Shroedinger equation then acquires the following form: 



i+> r (x,f) 
> r (x,t) 



[x,t) 

>,t) 



As-tfat) + V{x,t)s + ' i {x,t), 
As+'^z,*) + V{x,t)s-*{x,t), 
As + ' r (x,t) + V(x,t)s~' r (x,t), 
As-' r (x,t) + V(x,t)s + ' r (x,t). 



(14) 



We will enumerate the types (+,r),(+,i),(—,r),(—,i) by natural numbers 1,2,3,4 
respectively, and will apply to them arithmetic operations in Z/4Z everywhere including 
indices. We then have already no single density u, but the vector-column u, which compo- 
nents Uj j = 1, 2, 3, 4 are the densities of samples of types j. The rule for transformations 
of types must correspond to the turn around the center of coordinates in the definite 
direction, e.g., this rule must be cyclic: 1 — ► 2 — > 3 — ► 4 — > 1. The system (J14j) will 
be then equivalent to the following equation 



u = T(Au — qu), 



(15) 



where the matrix T, expressing our law of transformation of types, and the matrix g, 
inverting the signs of types have the form 



/ 



V 



o 
i 






1 \ 




0/ 



V 






1 





1 


0/ 



We now define the collective behavior giving the solution of equation ()15j) . Here it 
would not be sufficient to manage with only cyclic rule of transformation of types, because 
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in that case the matrix T would contain equal and nonzero diagonal elements, whereas in 
the matrix T all diagonal elements are zero. This is the reason of including tied photons 
to the collective behavior. 

We now have four samples of particle: 1,2,3,4 with the same dynamical characteristics 
and also the samples of tied photons. The samples of particle and photons are denoted by 
a J and 7 J correspondingly. Let each sample of charged particle of the type j = 1,2,3,4 
creates with the constant rate the samples of tied photons of the same type. These 
samples of tied photons move via diffusion and in the time frame St they transform to 
the samples of the same particle of the type j + 1. The type transformation then goes 
indirectly through tied photon samples. The main requirement is that the coefficient of 
diffusion for tied photon samples must be much nore than for the samples of particle: 
Pphot 3> Ppart- We can assume that the samples of particles do not move at all themselves. 
This allows to suppress the diagonal elements in the matrix of generalized diffusion that 
gives the needed mathix T. 

We then consider two areas Di and D 2 with the common border through which the 
diffusion goes. Here the photon samples diffusion goes much faster then the diffusion of 
samples of particle. In the change of density of any type of samples of particle in the tiem 
frame St two samples of this type make deposit: a) newly created as the transformation 
of type from the photon samples which either came from the neighboring areas, or were 
emitted in this area, and b) penetrated by the diffusion. By our supposition about 
intensities of diffusion, the deposit b) is negligible, and we can consider only deposit 
a). Applying the formula (Jll)) . we obtain that the income of samples of the type j in the 
unit time frame through the border is 

dQ = -pp ho t d _ dS, (16) 

where is the density of samples of parent type. We treat the potential V as the 
intensity of creation or annihilation of samples of any type. Using the reasoning from the 
derivation of the diffusion equation we come to the equation (JT3J). 

The introduced method of description of quantum evolution includes the classical 
dynamics of particles. It would suffice to have the number of photon samples A(|s r ' + — 
s r, ~| 2 + \s l,+ — s*' - ! 2 ) 1 / 2 , where A is constant which does not depend on the point of 
configuration space that prevent the useless storing in memory of the mutually canceling 
parts of the wave function. Let the carrier of wave function concentrate in some area 
-Do- The distribution of the phase as <p(x) = xp then lead to the same result as the 
mechanic movement of the swarm with the speed proportional to p along the vector p. 
Correspondingly, we can modify the rule for the transformation of photon samples to 
the samples of particle using the mean speed v av (x,t) of samples of particles which emit 
the samples of photon in the point x, t. The former main rule was that the photon 
samples of type j transform to the samples of particle of type j + 1 in time instant t + St, 
where St is constant. The new rule consists in that the photon samples transform to 
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the samples of particle of the types corresponding to the distribution of amplitudes of 
the form exp(5tip 2 / m) , where 5t is an arbitrary positive number in the segment (Q,8to), 
where the impulse p = v av m is computed basing on the mean speed of samples of particle 
which emit photon samples. It means that for the fixed St in this segment the fraction 
of the type j + k, k — 0, 1, 2, 3 equals to the corresponding amplitude in the exponential 
exp(5tip 2 1 m) . 

The new rule has no influence to the dynamics of swarm if the initial distribution of 
amplitudes has the form exp(ipx) and V = 0. In this case the movement is the uniform 
and has the constant impulse p. In the other cases the influence of this rule is determined 
by the value 5to: the less it is, the less the influence will be. If we fix the value 5to, we can 
find the influence of potential V. Let A^j(Ax) be the divergence of wave function values 
found by the old and new rules for the transformation of photon samples to the samples 
of particle in the point x + Ax. A^(Ax) will be then proportional to sin(grad V, Ax). 
The divergence is then maximal if the samples of photon move orthogonally to grad V. It 
substantiates the following agreement. We assume that a sample of tied photon transform 
to a sample of free photon if it propagates to the distance exceeding some limit Axq from 
the point of emission towards Ax, so that | sin(gradV, Ax)\ > e , e.g., in the direction 
close to the normal to the gradient of potential in which the considered particle moves. 
This agreement is not accurate because we have not defined the constants eo, Axo in 
the definition, but it can be reformulated in terms of QED if we introduce the vector of 
polarization of photon which is orthogonal to the photon impulse. Our agreement then 
will express in swarm terms the rule for finding of amplitude of fundamental process of a 
photon emission by charged particle (see (0)). 

4.2 Representation of Coulomb field 

The collective behavior defined below has one drawback: it does not explain completely 
the mechanism of photon exchange and thus it is not appropriate for the extrapolation 
to QED. Namely, the potential V(x, t) is used without explanation of its generation. It 
is sufficient for the creation of video film about many different particles interacting by 
Coulomb (or any other) law, but it cannot serve as a tool for description of emission 
absorption of the real photons. But the main deficiency is that the described method 
cannot be scalable due its huge computational complexity: for N samples to one particle 
and k particles one step of simulation requires of the order k 2 N 2 operations because 
with the assumption of instantaneous propagation of the field the interactions in all pairs 
of samples for the different particles must be considered. We show how to reduce this 
complexity to kN. For this we assume that the field propagates with the finite speed. In 
addition, our construction is true for the Coulomb field only. We introduce the mechanism 
of creation of Coulomb potential V(x, t) = \ X 2 XQ \ , by the center located in the point xq. 

Let the center emits samples g\,g2,... of new type free scalar photons, each of which 
is supplied by the sign of charge which has emitted it sign(gj), and the diffusion of these 
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samples go independently of the samples of charged particle; in the instant when they 
meet any sample of charged particle s a ' 5 , a G {+,—}, 5 G {r, i} they in turn emit 
the new sample of this particle of the type s a ' 6 , a' = —sign(s)sign(co)a. We assume 
that the free photon samples propagate much faster than the samples of tied photons. 
The propagation of the scalar photon density is thus determined only by the operator of 
kinetic energy E kin = — ^ and we obtain that with the constant scalar photon emission 
rate ^ in the point xq in the short time frame their density will be equal to Green function 
p{x) = — 47r 2m\ x -x \ °f ^ ne operator Ekin- Hence, after the transformation of scalar photon 
samples to the samples of probe particle located in the Coulomb field, we obtain the swarm 
behavior simulating one quantum particle in Coulomb field. The difference between free 
and tied photons is that the samples of free photons do not disappear, but create (emit) 
the samples of particles. 

The case of constant center x is idealized because it was assumed that the particle 
created the field has the infinite mass and thus does not reveal quantum properties. A 
real particle (for example, proton in the hydrogen atom) must be considered as a quantum 
particle as well, and, correspondingly must be represented by its swarm of samples d"' 6 . 
Each of these samples must emit scalar photon samples as above. A system of two charged 
particles is then represented by three swarms: 

{sj ,<5 , df, 7,<7,}, j = 1,2,...; (17) 

Standard quantum description of states in the form \E , (r 1 ,r 2 ) means that we consider 
the situation when \&(r 1 ,r 2 ) = ^/ + ' r (r 1 ,r 2 ) — ^~' r (ri,r 2 ) + i(^ +,l (r 1 ,r 2 ) — ^~ ,l (ri,r 2 ) 
^' 5 (ri, r 2 ) denotes the density of such pairs of samples of the form s a '' 5 ' , d a "' 5 " (relatively 
to the total number of such pairs over all a, 5), such that one sample has transformed 
from the scalar photon sample emitted by the other sample and (a, 5) = (a', 5') (a", 5") 
(product of pairs is defined by the natural way). The ordinary diffusion method Monte 
Carlo consists in that such pairs arise with intensity proportional to l/|ri — r 2 \. It is 
equivalent to that in the swarm approach any scalar photon sample g emitted by one 
particle samples in turn emit (or absorb, dependently on the charge signs) the sample of 
the second particle in the moment when g meets any sample of the second particle. The 
swarm representation of wave functions then corresponds to the computations along the 
Monte Carlo scheme in the framework of Hilbert formalism. 

We show how Coulomb interaction can be reduced to the swarm dynamics defined 
above. Let the swarm S corresponding to the charged particle be somehow divided to two 
parts: S' S", which we treat as the physically different quantum particles. Let at first 
these two swarms are separated by some spatial distance d. We consider the free photon 
samples which emit (create) the samples of particles and move themselves only by the 
diffusion mechanism. The density of these free photon samples will then be distributed 
by the law 1/r, where r is the distance passed by a photon sample from the point of 
emission, because 1/r is proportional to Green function of the diffusion operator A. If 
the photon samples emit (or absorb - dependently on the charge sign) the samples of the 
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other particle swarm when they meet its members, it is equivalent to the situation when 
they are initially distributed as 1/r, e.g., it is equivalent to Coulomb field. 

Hence the simple rule comes for finding of the approximation of Coulomb interaction. 
Coulomb interaction of two quantum particles is equivalent to the point wise interac- 
tion between their samples and the scalar photon samples emitted by the other particle 
samples, e.g., to the interaction which switches on when the coordinates of two samples 
corresponding to photon and particle coincide. We represent the quantum flow of one 
particle wave package through the tied photons transformations, which as the Coulomb 
free photon samples move via diffusion. Our representations of Coulomb field and free 
flow of wave package are then analogous. The accuracy of this approximation is limited 
by the supposition that the carriers of Coulomb field free scalar photons samples have 
much higher diffusion rate than the tied photon samples which determine the action of 
the diffusion operator with the imaginary coefficient: iA. The assumption of the point 
wise interaction makes swarm approach resembled to the density functional method, but 
this assumption in swarm representation will be valid in the three dimensional case only 
because only in this case Green function of the diffusion operator has the Coulomb form. 
In the other words, swarm method in the described form is applicable only to the Coulomb 
field in the three dimensional space. The supposition about point wise interaction radi- 
cally simplifies the practical realization of swarm approach for many particles as well as 
the representation of entangled states. Really, in any step of simulation the number of 
operations will be kN where k is the total number of particles, N - the total number of 
samples for one of them (and of photon samples emitted by its samples). The quadratic 
speedup comparatively with the search of all pairs "particle - particle" will be if we store 
in memory not the single photons but its density in all points of configuration space 
division. 

How to account the identity of particles in the swarm approach? The general idea 
is to treat the identical particles as the different areas of the same swarm of samples. 
These areas can be separated, say, by the sign of sample a so that for the positive sign 
we treat the sample as belonging to the swarm of the first particle, for the negative to 
the swarm of the second one. The Pauli principle is then guaranteed immediately. The 
more accurate consideration can be fulfilled if we introduce the spin for each sample of 
particles and make the intensity of interaction dependent on spins. 

4.3 Transformation of tied photons to free photons 

The swarm representation of ordinary quantum mechanics requires the introduction of tied 
photon samples which can be emitted by and can transform to the samples of charged 
particles. The description of QED effects requires the criterion for the transformation 
of tied photon samples to the free photon which can be detected. Such a criterion can 
be based on the visual representation of quantum dynamics only because it cannot be 
described inside QED itself. The following topological criterion can be used. We call 
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the division of swarm the situation when the spatial area occupied by this swarm is split 
to the different components of integrity: D\, D 2 , . . . , D^. One from the parts of swarm 
Si, S2, ■ ■ ■ , Sk, occupying these areas, say, S\ is then called the swarm of particle and all 
other parts form free photons of vector type (vector photons) about which we assume 
that they are emitted by this particle. The tied photons emitted by swarms S 2 , S 3 , . . . , Sk 
yet do not transform to the samples of charged particle but simply change their type 
correspondingly to the cyclic rule re, + — > im, + — > re, — — > im, — — > re, + 
and move with the constant speed corresponding to the impulse p of the charged particle so 
that the wave function found by the summing of photon samples amplitudes of four types 
approximately equals expiipx). Here we can separate the scalar photons propagating 
along the speed v of initial wave package and the vector photons propagated along the 
perpendicular to v. 

The transformation of free photon samples to the samples of charged particle goes 
along the rules dual to ones defined above. A free photon sample ph emits the sample of 
charged particle s when ph collides with a sample of s. This general receipt can be done 
in more details so that it will correspond to QED rules for the founding of amplitudes 
of emission and absorption of the virtual photons. For example, the probability of such 
transformation can be made proportional to the dot product of the polarization vector 
of ph and the impulse of the sample s it meets. The amplitude of emission (creation of 
photon sample by a sample of charged particle) and absorption (transformation of photon 
sample to a sample of charged particle) must coincide and must be proportional to the 
value 

ei(e,p), (18) 

where e and p are vector of photon polarization and impulse of charged particle, e is the 
charge of electron, i is imaginary unit. This rule of QED can be naturally reformulated for 
the swarm representation. This rule of QED can be naturally reformulated for the swarm 
representation. A photon sample arisen in the emission from some sample of charged 
particle of the type r must have the type (r + l)sign(e,p), where the intensity of this 
emission must be proportional to the module of this scalar product. This rule gives the 
swarm approximation of QED. The analogous procedure takes place for the particles with 
spin (see PJ). 

If the samples of charged particle are supplied with spins, we can introduce the spirality 
of free photon samples emitted by the samples s, of considered particle, so that the 
spirality of a photon sample ph uniquely determines the spin of sample s which emitted 
this photon sample. And vice versa, the photon spirality determines the spin of the sample 
of particle emitted by this photon sample when it meets any member of swarm s. The 
influence of the photon sample spirality will be opposite to the influence of spin of the 
sample emitting this photon sample. The closer is spin and spirality the easier will be the 
transformation from photon sample to the particle sample will be. We thus can simulate 
the effects connected with the interaction between spins of charged particles. 
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4.4 Representation of entangled states of many particle systems 



It was shown that the swarm method includes the dynamics of not entangled states. 
However, it is possible to extrapolate the swarm representation to the entangled states of 
the form (JBJ). Let we are given a system of k quantum particles. Its swarm description 
has the form: 

S\S 2 ,...,S k , 7, g, (19) 

where is the swarm corresponding to the particle j — 1, 2, . . . , k; 7, g are the common 
swarms of tied and free photon samples correspondingly. For each j we introduce the order 
on the swarm S 1 - 7 = {s{, s^, . . . , s 3 nj } and will consider some pairs of samples of the form 

sj, si as connected (z 7^ j, i,j = 1,2, ... ,k). We agree that if some sample disappears in 
the evolution then a sample of the other particle connected with it disappears as well. We 
then can simulate the entangled states of the form Let any state \&j be the coordinate 
position of the particle. The order on the swarms can be defined so that exactly those 
samples will be connected which correspond to the wave functions of the form ^ and ty™. 
We can also simulate the dynamics of such states by the changing of order on swarms. 

We can yet more simplify the representation if we require that exactly those pairs 
s\,s\ of samples are connected whose space-time coordinates coincide. At last we can 
suppose that the entanglement is transmitted by the photon sample from the sample 
which emitted it to the sample of the other particle. 

Coulomb interaction in the many particle system consists of point wise interaction 
between the samples of photons and particles. This interaction goes when the samples 
are located in the same point of space-time and with the intensity proportional to the 
probability of their collisions. Electromagnetic interaction goes through vector photons 
according the described rules. 

5 Complexity of the simulation of Coulomb interac- 
tion and the speed of its propagation 

The physical characteristic of field its speed of propagation is tightly connected with the 
complexity of simulation of the dynamics induced by this field. Let us consider once more 
swarm method of simulation of the dynamics of k particle system which representation 
has the form ()19|). In each step of the simulation we search all points of division of the 
configuration space, and for each point fulfill 

• Step of emission of the photon samples by the samples of particles located in this 
point, from the swarms Sj, 

• Step of the diffusion of the photon samples, located in this point, from the swarms 
7, g, and 
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• Step of emission of particle samples by the photon samples, located in this point. 

This method gives the complexity of simulation in each step of simulation of the order 
kN, where N is the number of all elements of configuration space. Let us compare this 
method with the direct method of simulation of the type each sample of i-th particle with 
each sample of j-th particle where the potential of interaction is proportional to l/r(i,j) 
in case of Coulomb field, and interaction propagates with the speed c. It requires of 
the order k 2 N 2 elementary operations to each step of iteration. To simulate one step of 
the direct method we need at most of the order c steps of swarm method. Hence, if c is 
finite, the complexity of collective behavior method has the order square root of the direct 
method complexity. But if c is so large that dr/c < dt, where dr, dt are the minimal 
sizes of the length and time in our model, our swarm method will work unsatisfactory. 
Really, if the field propagation in our representation is instantaneous, then the accuracy 
of our approximation will depend on the total number k, of particles, that violates the 
requirement of the scalability. We ignore the influence of particles located too far, but 
their summed influence in instantaneous interactions can be valuable. 

The requirement of the linear growth of computational recourses with the growth 
of k leaded to the limitation of the speed of field propagation in the representation of 
collective behavior. The finite value of c immediately gives the effect of delay. During the 
time required for the propagation of free photon samples, transmitting Coulomb field, the 
samples of tied photons change a little the shape of the particle wave function. It allows 
to simulate the relativistic effects in the method of collective behavior. 

We note that the large speed of free photon diffusion comparatively to the tied does 
not require the larger instantaneous speed of their moving. For this it would suffice to 
assume that that the intensity of creation of the samples of particles for the free photons 
is much smaller than for tied, but the intensity of their own creations is, in contrary, much 
higher. This relation makes possible to assume that the free photon samples transform 
to the samples of particles as the tied photon samples. 

The reasonings from above are valid for every field which propagation is determined by 
the local interactions of the type of cellular automata, as the diffusion method for Coulomb 
field. If algorithmic approach is universal and can embrace all the existed interactions 
in the common simulation algorithm, then we can roughly estimate the computational 
complexity of QED of electrons and atomic nuclei in comparison with the other inter- 
actions (strong, weak, and gravitational). If we suppose that the fraction of quantum 
non-locality in the common volume of computations is small, then the time of simulation 
r will approximately correspond to the physical time. The computational complexity of 
QED will be then proportional to the time needed for the processing of the free photon 
samples shift to the unit of length, relative to the time needed for the processing of the 
same shift of the samples of charged particles, e.g. the fraction of speeds of particles and 
photons: v par t/v p hot- The large value of the speed of light c in comparison with the typical 
speeds of charged particles serve as the indicator of the large computational complexity of 
the phenomena of not electrodynamic nature, which occupy the greater part of the time 
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of the universal simulating algorithm. 



6 Correspondence between swarm method and QED 

Swarm representation for many body systems with connected samples serves as a reduction 
of Hilbert formalism that is required by the algorithmic approach. This representation 
makes possible to build the films of processes which formulation in QED is incorrect: 
for example, the problem of capture of flying electron by a proton. Incorrectness of this 
problem in QED follows from that here the state of photon field is unknown, e.g., there are 
no photons at all. Hence, this problem cannot be solved through the scattering matrix. 
At the same time the inverse process - photo effect has the standard description in QED. 
Factually, the reason of incorrectness of such problems in QED is the general, and it 
follows from the principal impossibility to glue quantum mechanics with the classical. 
The difficulty is not in that it is impossible to join these two techniques by the unit set of 
formulas it would join easily in the framework of the common simulating algorithm. The 
difficulty is that there is no criterion which determines for a given system whet technique 
to apply for it: classical or quantum. It explains the importance of the wave function 
collapse problem. The creation of a scalable quantum computer would not solve this 
difficulty but only reduces its sharpness because a quantum computer makes possible to 
extrapolate the existing methods to the wider class of problem without their substantial 
modification. But the absence of clear progress in the technologies of quantum processors 
gives grounds to suppose that that this difficulty is solved in Nature by some mechanism 
which immediate actions we observe in experiments but the explanation of it does not 
require to address to extremely complex objects 6 . It gives the new sense to the attempts 
to find the modifications of QED which would be free of this difficulty, particularly, the 
modification based on the collective behavior. 

This task is substantially harder than the reduction of ordinary quantum mechanics 
to the collective behavior. At first, here we deal not with the logically simple and reliable 
quantum mechanics of Shoedinger equation, but with the logically incomplete quantum 
electrodynamics where such things as divergence of rows take place. At second, in the 
formalism we speculate on the incorrect in QED problems must have the certain solution. 

Factually, the passage "quantum mechanics of many bodies j 1 classical mechanics" 

must become the correct and monosemantic, that augurs much radical logical break with 
Hilbert formalism than the passage to the collective behavior in ordinary quantum me- 
chanics. Classical algorithms, which earlier served at most as the replacement of the 
calculations by formulas, will now come in the conceptual part of description, and it is 
unavoidable if we intend to reach our main goal. It worsen dramatically perspectives to 
find the correspondence with QED. If for the ordinary quantum mechanics such corre- 

6 for example, to the consciousness of an observer (compare with 6 ). The algorithmic approach, 
developed in this work is the search of such mechanism. 
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spondence can be lead to the form of mathematical theorems 7 , then for the many body 
case with the field even the search of formulations of such theorems represents problems. 
It is impossible to substantiate completely our model in the framework of QED, because 
the divergence with QED is laid in the starting point of swarm model. The reasonable 
form of checking of all versions of algorithmic formalism including the method of collective 
behavior would be the description of the concrete processes and the comparison with the 
standard when the both descriptions are possible. The algorithmic approach pretends in 
perspective to the including to the physics the scientific areas which traditionally belong 
to the other sciences (chemistry), and the algorithmic representation of the dynamics 
presumes the visualization. But the structure of simulating algorithms must be simple 
enough to allow the comparison with the standard quantum description when the last 
is possible. It touches the scattering problems that are computed with high accuracy 
including relativistic corrections, etc. The requirement of good approximation in these 
cases is necessary for the development of our approach 8 . 

7 Conclusion 

We have defined the algorithm for the simulation of many particle quantum dynamics 
accounting the simplest EPR type of entanglement. This algorithm in the elementary 
form includes the Coulomb interaction only, but it can be supplied with the exchange 
of vector photons leading to the effects of QED. The main advantage of this algorithm 
is its scalability when the total number of particles grows. The further development of 
this algorithm is connected with the visual representation of the many particle dynamics. 
The visualization is the principal tool of the algorithmic approach which makes possible 
to generalize the method of theoretical physics on the complex many particle systems 
studied in chemistry and biochemistry. 

7 Feynman path integrals give the equivalent form of quantum mechanics which does not touch its 
conceptual basement. But in the exact proof of this fact some suppositions must be introduced that are 
not evident beforehand: for example, about the uniform distribution of linear parts of paths and the 
complete substantiation requires certain efforts. But in case of new formalism which does not completely 
agree with Hilbert these difficulties grow many times. For the ordinary quantum mechanics of one particle 
this is not felt at all. 

8 The method of collective behavior cannot compete successfully with the analytical methods in any 
cases, because the real total number of samples can be much more than can be stored in the memory of 
our computers. But we expect that this method can substantially widen the borders of processors which 
we can simulate somehow. However, sometimes the methods of swarm type can beat analytical methods 
in the precision. The good example is the diffusion Monte Carlo method for the ground states of two 
electron system (for example, atom of helium or molecule of hydrogen). 



19 



References 

[1] [FejR.Feynman, The theory of fundamental processes, Caltech, NY, 1961. 

[2] [FH]R.Feynman, A.Hibbs, Quantum mechanics and path integrals, McGraw-Hill 
book company, NY, 1965. 

[3] [OzjY.Ozhigov, Algorithmic approach to quantum physics, lanl e-print 
quant-ph/0412196, journal version in Mikroelektronika, 2006, vol. 36, 1, pp. 
61-77. 

[4] [SejI.Semenihin, Y.Ozhigov, Algorithmic approach to quantum theory 2: Amplitude 
quanta in diffusion Monte Carlo method, Proceedings of SPIE, 2006, 6264. 

[5] [KB]P.Kurakin, H. Bloom, G.Malinetskii, Conversational (dialogue) model of quan- 
tum transitions, lanl e-print |quant-ph/ 0504088 , Proceedings of SPIE, 2006, 6264. 

[6] [MejM.Mensky, Quantum measurement: decoherence and consciousness, UFN, 2001, 
vol. 171, 4, p. 459. 

[7] [Vl]Vladimirov, Equations of mathematical physics (rus.), Moscow, Nauka, 1974. 



20 



